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ABSTRACT 

Between 7 March 2002 and 15 June 2002, intensive X-ray observations were 
carried out on the extreme BL Lac object H1426-I-428 with instruments on board 
the Rossi X-ray Timing Explorer (RXTE). These instruments provide measure- 
ments of H1426-I-428 in the crucial energy range that characterizes the first peak 
of its spectral energy distribution. This peak, which is almost certainly due to 
synchrotron emission, has previously been inferred to be in excess of 100 keV. By 
taking frequent observations over a four-month campaign, which included ~450 
ksec of RXTE time, studies of flux and spectral variability on multiple timescales 
were performed, along with studies of spectral hysteresis. The 3-24 keV X-ray 
flux and spectra exhibited significant variability, implying variability in the lo- 
cation of the first peak of the spectral energy distribution. Hysteresis patterns 
were observed, and their characteristics have been discussed within the context 
of emission models. 

Subject headings: AGN, BL Lac, H1426-h428, X-ray 

1. Introduction 

Blazars, including flat spectrum radio quasars (FSRQs) and BL Lacs, represent the 
most extreme examples of active galactic nuclei (AGN) due to their rapid variability and 
the extreme energies to which their spectra extend. One of the most exciting developments 
in extragalactic astrophysics during the last decade has been the detection of more than 65 
AGN in the energy range from 30 MeV - 30 GeV (Hartman et al. 1999) by the EGRET 
experiment on the Compton Gamma Ray Observatory (CGRO). The detection of some AGN 
(e.g. Mrk 421, Mrk 501, lES 2344+514, PKS 2155-304 and most recently H1426+428 and 
1ES1959+650) above 300 GeV by ground-based air cherenkov telescopes (ACTs) has rein- 
vigorated the study of these objects (Catanese and Weekes 1999; Horan et al. 2002; Holder 
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et al. 2003). The dominant radiation from blazars is widely believed to arise from rela- 
tivistic jets that are viewed at small angles to their axes (Blandford and Koenigl 1979). 
EGRET and Whipple observations of high fluxes above 1 GeV and 300 GcV (respectively), 
the short-term variability of the gamma-ray emission (Kniffen et al. 1993; Mattox et al. 
1993; Quinn et al. 1996; Buckley et al. 1996), and constraints on gamma-ray opacity due to 
pair-production with low energy photons are all pieces of evidence leading to the conclusion 
that the gamma-ray emission is also beamed and that jets are involved. 

The broadband radiation spectrum of blazars consists of two parts: a synchrotron spec- 
trum that spans radio to optical-ultraviolet wavelengths (and to X-rays for high-peaked 
objects) and a high-energy part that can extend from X-rays to gamma rays. Observation- 
ally, the spectra appear to have two distinct bumps when plotted in a i/fi, representation 
(Fossati et al. 1998; Ghisellini et al. 1998). The emission in the first peak of the spectral 
energy distribution (SED) is generally considered to be synchrotron emission from relativis- 
tic electrons. The most widely invoked models which attempt to explain the higher energy 
emission of the second SED peak fall into the category of leptonic models. These leptonic 
models posit that the X-ray to gamma ray emission is produced by inverse Compton scat- 
tering of lower energy photons by relativistic electrons within a narrow jet. The low energy 
photon fields could arise from synchrotron continuum photons within the jet (e.g. Koenigl 
(1981)), or they could arise from ambient photons from the accretion disk which enter the jet 
directly (e.g. Dermer, Schhckeiser, & Mastichiadis (1992)) or after scattering or reprocess- 
ing (e.g. Sikora, Begelman, & Rees (1994)). In addition to these leptonic models, hadronic 
models also attempt to explain the second bump of the SED. These models involve proton- 
initiated cascades (e.g. Mannheim (1993)) and/or proton synchrotron radiation (Muecke 
and Protheroe 2001; Aharonian 2000), as well as synchrotron emission from secondary 
muons and pions(Muecke et al. 2003). 

The gamma-ray and simultaneous X-ray observations of blazars have placed strong con- 
straints on these models (von Montigny et al. 1995; Macomb et al. 1995; Buckley et 
al. 1996) by constraining the properties of the emission regions (size, Doppler boosting 
factor, magnetic field, electron cooling time, total injection energy, ...) and by comparing 
fiaring timescales at different wavelengths; but so far the observations cannot definitively 
reject either hadronic or leptonic models of the high energy emission. X-ray observations 
are particularly important for observations of high-peaked BL Lacs since the peak of the 
synchrotron component lies in the band covered by X-ray detectors (e.g. Beppo-SAX ob- 
servations of Mrk 501 (Plan et al. 1998) and H1426+428 (Costamante et al. 2001)). The 
determination of the X-ray peak is highly constraining to the suite of models for the emission 
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from these objects, as are the timescales of the X-ray variabihty. 

Further constraints to emission mechanisms may be obtained by exploring the charac- 
teristics of spectral hardness-intensity diagrams (HIDs) in the X-ray energy band, with the 
purpose of searching for hysteresis patterns. The presence of hysteresis loops and the handed- 
ness of such loops can be used to constrain the cooling and acceleration/injection timescales 
of the relativistic electron population (Bottchcr and Chiang 2002; Li and Kusunose 2000; 
Kirk, Ricgcr, and Mastichiadis 1998). Spectral hysteresis has been observed for several 
high-pcakcd BL Lac objects, including Mrk 421 and PKS 2155-304 (Takahashi et al. 1996; 
Kataoka et al. 2000; Zhang 2002). Another set of observations of Mrk 421 using XMM 
reported no significant hysteresis (Sembay et al. 2002), however this result appears to be 
inconclusive due to the lack of coverage throughout the entire cycle of flux rise and subse- 
quent decay coupled with the low amplitude variation of the one complete rise and decay. 

H1426-I-428 is an X-ray selected high-peaked BL Lac with a relatively strong X-ray 
flux. It was flrst characterized as a BL Lac, due to its lack of prominent optical emission 
hues, by Remillard et al. (1989). Based upon observations by Beppo-SAX up to 100 keV, 
H1426-I-428 is an example of an extreme high energy peaked BL Lac (HBL) with the peak 
of the synchrotron emission lying above 100 keV (Costamantc ct al. 2001) during that 
observing campaign. This fact led Costamante et al. (2001) to speculate that H1426-I-428 
is a prime candidate for TeV gamma-ray emission; a speculation recently confirmed by the 
Whipple collaboration (Horan et al. 2002; Petry et al. 2002). Characterizing the X-ray 
properties of this extreme object is important, especially due to the presence of this TeV 
emission in spite of a large redshift (z=0.129). Previous X-ray observations presented by 
Sambruna et al. (1997) and Madejski et al. (1992) using ASCA, ROSAT, and BBXRT 
(at different epochs separated by as much as 15 years) have shown evidence for spectral 
variability from one observation to the next, as well as evidence for spectral curvature and 
a spectral absorbtion feature at ~0.6 keV. The results presented here represent a fairly 
continuous and deep exposure obtained for the purpose of exploring the details of the X-ray 
variability on multiple timescales. 

2. Observations 

X-ray observations were carried out between 7 March and 15 June of 2002. These 
observations were part of a multiwavelength campaign (Falcone et al. 2003), including 
radio, optical. X-ray, and gamma-ray energy bands. This report will describe the X-ray 
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observations (multiwavelength observations will be reported on in a separate publication). 
X-ray observations were performed by PCA and HEXTE on the RXTE satellite. The results 
reported here are based on PCA data. The HEXTE data, which can only be analyzed 
on significantly larger timescalcs due to sensitivity constraints, will be described in the 
forthcoming multiwavelength paper. Typical RXTE exposures were between 1-4 hours per 
night for 3 weeks bracketing the new moon of each month during the 4 month long campaign. 
This resulted in ~450 ksec of RXTE observations. 



2.1. Instrument 

The PCA instrument is composed of five proportional counter units (PCUs) and has 
a total collection area of 6500 cm^ when all counters arc operating. However, during the 
observations discussed in this paper, only 2-3 of the proportional counters were on during 
normal operation. For nearly all of this campaign, the proportional counters known as PCUO 
and PCU2 were collecting data, while the other PCUs were dormant. For consistency during 
data reduction, a cut was placed on the data to select only the PCUO and PCU2 data. Data 
from these PCUs were extracted independently of one another and they were combined by 
jointly modelling the spectra at the time of analysis. The best statistics were obtained by 
using the 1st Xenon layer of each PCU and excluding layers two and three. Cuts were also 
made to eliminate data that included effects due to passage through the South Atlantic 
Anomaly, as well as data for which the electron rate was anomalously high. An elevation 
cut (>10deg) was also imposed. 

The extraction and processing of the data was done using ftools 5.2. A response matrix 
was calculated, using the latest available software (PCARSP V8.0), individually for each 
PCU with the same configuration used for the data extraction process. The weak source 
background model provided by RXTE/GOF was used for background subtraction. This 
model appeared stable and reproduced the observed spectrum at high energies, where the 
spectrum is expected to be background dominated. It is worth noting that the background 
model worked well for PCUO, which no longer has an operating propane veto layer. The 
model background for PCUO converged to the observations at high energies and the rate 
variations in PCUO mimicked those of PCU2 which has an operating propane veto layer. 
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2.2. Analysis 

After extraction of the data, XSPEC 11.2.0 was used to analyze and apply models to 
the extracted spectra. The nominal energy range of the PCA is 2 to ~60 keV, but since 
the response of the instrument is poorly characterized in the lowest energy bins, only data 
above 2.9 keV is used. Due to the combination of faUing source spectra for H1426+42.8 
and sensitivity as energy increases, there is no significant measurement above 24 keV. Based 
upon these two bounds, all analysis is confined to the region between 2.9 and 24 keV, and 
all other energy bins are excluded. The 3.8-5.4 keV energy region is also excluded from the 
analysis since there is a known Xenon absorbtion feature within this energy region that is 
not modelled well by the instrument response. The spectra are then fit to data that are 
binned into 1— day bins, as described below, and the fits are evaluated using a analysis 
and an inspection of the residuals. Count rates in particular energy regions are obtained 
prior to application of spectral models in order to evaluate variability in the lightcurve. Flux 
is derived based upon the parameters of the model that fits the spectral data. 



3. Results 

3.1. Flux and Flux Variability 

A lightcurve showing the X-ray rate as a function of time in different X-ray energy bands 
is shown in figure 1. This data has been binned into one day time intervals beginning with 
the first day of the campaign, and all error bars are la statistical errors. The 2.9-10 keV and 
the 10-24 keV energy bands have been evaluated independently of one another and plotted, 
along with the entire 2.9-24 keV band. The 2.9-24 keV flux observed by PCA varied between 
about 1.3x10"^^ to 6.6x10"^^ erg cm~^ s~^, and the mean 2.9-24 keV flux was 3.6x10"^^ erg 

— 2 —1 

cm s . 

Variability is evident in all of the three bands. Based on this analysis, which uses bins of 
1-day duration, variability studies are limited to the timescale range bracketed by 1-day 
and 4-months. There appears to be long term variability as the flux changes throughout 
the entire campaign with timescales in excess of two weeks, and there also appears to be 
variability on timescales as short as 1-2 days. While the variability was signiflcant, especially 
on long timescales, there were no particularly large flares relative to those observed for other 
TeV emitting BL Lacs, such as Mrk 421, Mrk 501, and lES 1959+650. 
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Fig. 1. — X-ray count rate per PCU as a function of time during Mar-Jun of 2002 in three 
energy bands. 
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3.2. Spectra and Spectral Variability 



For each set of data in a one day time interval, a model of the spectra was derived. It 
was found that a power law with galactic absorbtion, column density Nh = 1-4 x 10^'' (Stark 
et al. 1992), led to an excellent fit based upon a x^/dof analysis. (In fact, rather than 
hovering at ~1, the jdof is actually consistently less than one, due to the fact that the 
background rate error bars, which are added in quadrature to the source rate error bars, are 
conservatively large. This occurs since the error bars on the background are assigned based 
on the counts in the observation, rather than the counts used to construct the background 
model, which is a sum of many observations with far more numerous counts.) The power 
law spectra are of the form: 



where the units of the normalization constant C are photons keV'^cm'^s'^, cr{E) is the 
photo-electric absorbtion cross section given by Morrison and McCammon (1983), and F 
is the power law spectral index. Allowing the absorbtion column density to vary as a free 
parameter did not significantly effect the fit so it was fixed to the galactic value quoted 
above. An example of a fit using this data is shown in figure 2. 
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Fig. 2. — An example of a model fit to a typical day in the middle of the campaign (6 May 
2002), using the power law with galactic absorbtion model to jointly fit PCUO and PCU2 
data. Residuals of the fit (in terms of sigma) are shown in the bottom panel. 
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The power law spectral index derived using the above procedure is shown as a function 
of time in figure 3, along with /dof for the particular fit. All error bars are la. The 
error bars for the 2.9-24 keV flux were calculated by applying simple error propagation to 
the la error bounds produced by XSpec for the power law spectral index and normalization 
constant. It can be seen from this plot that there was significant spectral variability during 
the campaign. The spectral index varies over the range from 1.46±0.05 to 2.03±0.09. The 
X-ray flux in the 2.9-24 keV band has also been derived based upon the spectra that were 
fit to the data. This flux is shown in the top panel of figure 3. 
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Fig. 3. — Top panel: Xray flux as a function of time in the 2.9-24 keV band; bottom two 
panels: X-ray power law spectral index as a function of time along with the associated 
X^/dof for the fit. Error bars are la. 



3.3. Variability Correlations and Hysteresis 



Over the time frame of the entire campaign, there is no tight correlation between the X- 
ray flux and the power law spectral index. The correlation coefficient between the 2.9-24 keV 
flux and the spectral index has been calculated using the 1-day timescale data throughout 
the entire four-month campaign, and its value is 0.51. By simple inspection of the plots 
in figure 3, one can see that the spectral variability does not have a consistent relationship 
with the 2.9-24 keV fiux variabifity over long timescales. A plot of the spectral index vs. 
fiux is shown in figure 4. While, on average throughout the campaign, the flux increases 
as the spectrum becomes harder, there is a great deal of scatter on the HID for the entire 
campaign. This plot includes flux variations averaged over long timescales so it is certainly 
not reflective of the source characteristics during a single flare event. This lack of strong 
correlation does not preclude a relationship during individual flare timescales. 
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Fig. 4. — Hardness-intensity diagram for the complete four-month long data set. A hnear 
fit to the data has been overlayed on the figure to guide the eye. 
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Table 1: Definitions of time regions during periods of fiux variation 



Region 


Time 


Characteristics 


A 


1-12 May 


long rise, fast decay 


B 


12-16 May 


2-day rise, 2-day decay 


C 


30May-5Jun 


4-day rise, 2-day partial decay 


D 


5Jun-15Jun 


long rise, 2-day partial decay 


E 


7-11 Mar 


decay phase only (no observations during rise) 


F 


3-20 Apr 


long rise (no observations during decay) 


G 


30May-15Jun 


long rise (no observations during most of decay) 
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Hardness-intensity diagrams have also been produced from the data during individual 
periods of flux variation in order to search for evidence of hysteresis effects, which could 
be used to evaluate the cooling and acceleration timescales. Based upon the behavior of 
the flux, several time periods have been evaluated independently, as defined in table 1 and 
shown in figures 5 and 6. The time periods evaluated in figure 5 are those that have a rise 
and a subsequent decay in flux observed by RXTE-PCA. For three of these time periods, a 
well-defined spectral hysteresis is observed with the trace following a clockwise pattern in 
the representation plotted. For the fourth (region B), there is also a consistent hardness-fiux 
relationship observed, which tends to harden the spectrum with increasing fiux, but the 
handedness of the HID can not be determined within the resolution of the measurements. 
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Fig. 5. — Hardness-intensity diagrams for regions A, B, C, and D (clockwise from top-left) 
as defined in table 1, which all exhibit a rise and subsequent decay of flux. Arrows have 
been inserted to show the direction of temporal evolution. 
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The time periods evaluated in figure 6 are those which do not have a well-defined rise 
and decay, but do exhibit variable flux. Region E exhibits a fast decay, but there were no 
observations during the rise of the flare. The HID exhibits a softening spectra with the 
decrease in flux. Region F is characterized by a long rise in flux over a time period of >2 
weeks, but there are no observations during the decay of flux. This HID diagram exhibits 
an overall hardening of the spectrum as the flux increases, but there is some scatter in the 
data points. Hysteresis properties cannot be determined due to the lack of data during the 
presumed decay of the flux. Region G is similar to region F in that it is a long rise of fiux 
over a period of ~2 weeks, but it does turn over and begin to decay in the last two days. 
This HID exhibits spectral hardening with increased flux. An overall clockwise hysteresis 
appears to be evident for region G, however this is not well determined due to two factors; the 
complete long-timescale flux variation is not observed throughout the entire decay, and the 
short-timescale flux variations make it difficult to evaluate the characteristics of the longer 
timescale variability. The two smaller HID loops associated with regions C and D flaring 
can also be seen in this HID. The effect of relatively short term variability combined with 
long-term variability is evident in the HID diagram for region G. 
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Fig. 6. — Hardness-intensity diagrams for regions E, F, and G (left to right) as defined in 
table 1, which only exhibit a fraction of a flux rise/decay cycle. Arrows have been inserted 
to show the direction of temporal evolution. 



-18- 



4. Discussion and Conclusions 

This long campaign was designed to continuously monitor H1426+428 at X-ray energies 
(and at other wavelengths that will be reported on in a forthcoming paper) in an effort to 
explore variability at multiple timescales. In spite of the fact that there were no large flares, 
relative to those frequently observed for other TeV emitting AGN, the campaign successfully 
observed variability with doubling timescales ranging from ~1 day to >2 weeks. In several 
cases, the short timescale flaring was observed on top of the long timescale gradual varia- 
tions. This type of variability complicates studies of emission mechanisms and should be 
considered in general, especially for studies that include non-contemporaneous multiwave- 
length observations. 

Throughout the entire campaign, a simple power law with galactic absorbtion was able 
to flt the data well between 2.9 and 24 keV. Prior observations performed by Sambruna et al. 
(1997) between 1985 and 1994 with ASCA and BBXRT led to similar spectra at some times, 
but at other times there were some notable differences. At some times, the spectra from 
Sambruna et al. (1997) could be flt by a single power law, but at other times a broken power 
law that became softer at higher energies was inferred. For the part of the spectrum between 
the break energy and ~10 keV, the spectra had a typical power law index of ~2.3 above a 
break energy of ~2 keV for the observations described by Sambruna et al. (1997). This is 
somewhat softer than any of the spectra derived during the present campaign. The observa- 
tions reported by Sambruna et al. (1997) all show 2-10 keV fluxes that fall within the range 
of fluxes observed during this campaign, whereas the spectra of these historical observations 
are softer than the more recent observations reported in this paper and in Costamante et 
al. (2001). This is most easily interpreted as time variability, rather than interpreting it 
as a contradiction between the two results, when consideration is given to the BeppoSAX 
observations of Costamante et al. (2001) which also leads to a different spectral index and 
when consideration is given to the observed variability during this four month campaign. It 
is also important to note that this has involved a comparison of different instruments with 
different band passes. Another interesting feature of the BBXRT observations reported by 
Sambruna et al. (1997) was the evidence for a spectral line at ^0.6 keV, which would imply 
the presence of absorbing material. No spectral lines were observed in the 2.9 to 24 keV 
band covered by the observations presented in this paper (PCA can not observe down to 0.6 
keV where the spectral feature was previously observed). 

The power law spectral index was observed to vary throughout the range bounded by 
1.46±0.05 and 2.03±0.09. If one interprets this as being due to a shift in the location of 
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the first peak of the spectral energy distribution, then it imphes that the peak is sometimes 
in excess of 100 keV, while at other times it falls into the 2.9-24 keV region observed dur- 
ing this campaign. While past observations of H1426+428 have implied a high first peak 
in the spectral energy distribution (SED) (Costamantc ct al. 2001), there have also been 
observations from which a spectrum with r>2 was derived (Sambruna et al. 1997). Based 
on this campaign it can be seen that significant spectral variability is present on multiple 
timescales. The average spectral index was generally harder than a fiat i/f,^ spectrum, which 
imphes that the first peak of the average SED was probably in excess of ~100 keV (PCA 
alone cannot derive the actual peak location for these cases, but past observations (Costa- 
mante et al. 2001) with BeppoSAX have shown the peak to be in excess of ~100 keV 
when the spectral index is ~1.9 at PCA energies). At other times the spectral index was 
fiat (r ~2), which implies that the first peak of the SED was in the observed 2.9-24 keV 
region. It is interesting to note that this movement in the first peak of the SED during this 
campaign is not accompanied by any measurable fiuctuations in the TeV emission from this 
known TeV source (Falcone ct al. 2003), but the probability of detecting such associated 
TeV fluctuations from such a weak source has yet to be evaluated. The overall spectral 
energy distribution and its implications will be discussed in more detail in a forthcoming 
article that will incorporate radio, optical. X-ray, and VHE gamma-ray observations carried 
out contemporaneously during this campaign. 

The longest flux variability timescale observed was the long rise and short decay in 
June (region G) . In spite of the fact that only a few data points were available for the decay 
phase, some hysteresis was evident. The clockwise loop approached closure, as shown in the 
third panel of figure 6. Although it is only marginally significant, it is interesting to note 
the small loop in the middle of the HID diagram, which coincides with the shorter timescale 
fiux variability that occurred in the midst of the long rise in June. This behavior illustrates 
the need for sufficient sampling over complete fiaring timescales. It is probably this effect of 
hysteresis within hysteresis due to multiple flares, as well as flares that are non-simultaneous, 
that causes the scatter on the overall HID plotted in flgure 4. In spite of this scatter, it has 
been shown that there is a tighter relationship between flux and spectral hardness when each 
flux variability episode is studied independently on its proper timescale. 

Spectral hysteresis has been observed during several periods of flux variability during 
this campaign. For an unambiguous observation of spectral hysteresis, it is necessary to 
observe during a large fraction of both the decay and the rise of a flare. For three of the 
four cases in which this condition was met during this campaign, spectral hysteresis loops 
tracing a clockwise pattern were observed (given the axis orientation shown in flgure 5). In 
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the fourth case (region B), no hysteresis and associated orientation could be observed, but 
the spectrum did harden as the flux increased during the rise and decay of flux. In the 
past, two other high-peaked BL Lac AGN that are known TeV emitters (Mrk 421 and PKS 
2155-304) have exhibited similar clockwise loops in their hysteresis plots (Takahashi et al. 
1996; Kataoka et al. 2000). As mentioned previously, Sembay et al. (2002) did not observe 
hysteresis with XMM observations of Mrk 421, however these observations did not sample 
the rise and fall of the flux variability so it is reasonable to expect a null result in that 
case. The one flare that was completely observed by Sembay et al. (2002) did display the 
characteristic spectral softening with decreasing flux, but no loop could be discerned from 
this low amplitude variation. Another observation of Mrk 421 using BeppoSAX (Zhang 
2002) showed evidence of spectral hysteresis with the orientation in the counter-clockwise 
direction. So, for at least one object, the hysteresis characteristics can change from flare 
to flare. Since the hysteresis probes the timescale of the cooling processes as a function 
of energy, it appears as though some flares have a cooling time that acts faster at higher 
energies, consistent with synchrotron cooling, while other flares exhibit cooling times that are 
longer. For the H1426-I-428 observations presented here, the generally clockwise hysteresis 
curves are consistent with the scenario in which the X-ray flux is dominated by synchrotron 
emission since the synchrotron cooling timescale is shorter at higher energies. In contrast, 
if the hysteresis loops were oriented in the opposite direction, one would expect the source 
to be in a regime in which the cooling and acceleration times were approximately equal, 
according to models of Kirk, Rieger, and Mastichiadis (1998). 

By performing a study of the X-ray hysteresis, one can also explore the potential of 
various acceleration mechanisms to produce the TeV emission. Typically, clockwise hystere- 
sis patterns (using the axis orientation shown in preceding flgures) are interpreted as a soft 
time lag in X-ray emission, while counter-clockwise orientation is interpreted as a hard lag 
(Kirk, Rieger, and Mastichiadis 1998; Kataoka 2000). A soft lag at X-ray energies, which 
is interpreted as the energy dependence of the synchrotron cooling time, is characteristic of 
synchrotron self Compton models (Takahashi et al. 1996). A hard lag is not consistent 
with some standard one-zone synchrotron self Compton models, while models such as that 
of Kirk, Rieger, and Mastichiadis (1998), which incorporates a time-dependent propagating 
shock that accelerates electrons in the region of the shock front and allows for the case where 
Tcooiing ~ Taccel, Can bc Consistent with a hard lag. The observations presented here are gen- 
erally in agreement with standard models that predict a soft lag, but the lack of hysteresis 
from the flare in region B does not favor simple one-component models during that partic- 
ular time period. However, the model of Li and Kusunose (2000) can reproduce hysteresis 
loops that are not clockwise (and presumably intermediate situations with no obvious loop 
orientation) by increasing the injection energy to the point where synchrotron self Compton 
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losses begin to dominate the electron cooling. The TeV emission during this entire campaign 
was not in a high flaring state so firm conclusions regarding TeV acceleration during any 
short time period cannot be made. 

Detailed modeUing can be used to interpret the observed hysteresis patterns. Bottcher 
and Chiang (2002) have worked on such modelling of low-peaked BL Lacs, and they have 
found that external Compton and synchrotron self Compton mechanisms will produce dif- 
ferent signatures at X-ray energies. However, this study was aimed at low-peaked BL Lacs, 
rather than high-peaked BL Lacs. We have found that a comparison to these models is 
inconclusive, probably due to the significant difference in the parameter space sampled by 
X-ray observations since low peaked objects are viewed in an energy region much higher 
with respect to the synchrotron peak, relative to high-peaked objects. Li and Kusunose 
(2000) did model the hysteresis characteristics of high peaked BL Lacs, but they confined 
themselves to only synchrotron self Compton mechanisms. They did find that the hystere- 
sis characteristics are sensitive to total injection energy leading them to conclude that it 
would be difficult to draw any conclusions based upon hysteresis observations. However, 
for high-peaked BL Lac observations at energies at or below the synchrotron peak (such 
as the observations presented here), Li and Kusunose (2000) found that clockwise orien- 
tation could be expected when the injection energy was low and synchrotron losses were 
dominant, and the opposite orientation could be expected when the injection energy was 
high and inverse Compton emission was dominant. This supports the conclusion that for at 
least some of the fiares observed during this campaign (those with clockwise orientation) the 
synchrotron losses were the dominant cooling mechanism. More detailed modelling, which is 
beyond the scope of the current paper, of X-ray spectral characteristics that can be used to 
differentiate between various high energy emission models for high-peaked BL Lacs is needed. 
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